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Abstract. 



In this paper we try to bridge breakthroughs in quantitative sociology /econometrics pioneered during the last 
decades by Mac Fadden, Brock-Durlauf, Granovetter and Watts-Strogats through introducing a minimal model 
able to reproduce essentially all the features of social behavior highlighted by these authors. 
Our model relies on a pairwise Hamiltonian for decision maker interactions which naturally extends the multi- 
populations approaches by shifting and biasing the pattern definitions of an Hopfield model of neural networks. 
Once introduced, the model is investigated trough graph theory (to recover Granovetter and Watts-Strogats 
results) and statistical mechanics (to recover Mac- Fadden and Brock-Durlauf results). Due to internal symmetries 
of our model, the latter is obtained as the relaxation of a proper Markov process, allowing even to study its out 
of equilibrium properties. 

The method used to solve its equilibrium is an adaptation of the Hamilton- Jacobi technique recently introduced 
by Guerra in the spin glass scenario and the picture obtained is the following: just by assuming that the larger 
the amount of similarities among decision makers, the stronger their relative influence, this is enough to explain 
both the different role of strong and weak ties in the social network as well as its small world properties. As a 
result, imitative interaction strengths seem essentially a robust request (enough to break the gauge symmetry in 
the couplings), furthermore, this naturally leads to a discrete choice modelization when dealing with the external 
influences and to imitative behavior a la Curie- Weiss as the one introduced by Brock and Durlauf. 



^e-mail:adriano .barraOromal . inf n. it 
^e-mail:elena. agliariOf is .unipr . it 



1 Summarizing some main results of quantitative sociology 



In recent years there has been an increasing awareness towards the problem of finding a quantitative 
way to study the role played by human interactions in shaping behavior observed at a population level, 
ranging from the context of pure sociology to the one belonging to economic sciences. The conclusion 
reached by all these studies is that mathematical models have the potential of describing several features 
of social behavior, among which, for example, the sudden shifts often observed in society's aggregate 
behavior [30] [3S] [13] , and that these are unavoidably linked to the way individual people influence each 
other when deciding how to behave (the phase transitions in the language of thermodynamics |18]). the 
whole suggesting a promising potential application of disordered statistical mechanics to this field of 
research I2][l2[12j. 

Here we summarize what we understood as real breakthroughs in these analysis, highlighting two main 
aspects dealing with topological investigations on the structure of the graph built by social interactions 
and the kind of interactions themselves. Namely, the discovery of the fundamental role of weak ties in 
bridging different communities (due to Granovctter [T31 [Ml [13 [IS]) and the "small world" feature of 
the social structure (obtained by Watts and Strogatz (TS] [38] [39] ) for the first analysis and the discrete 
choice of decision makers in econometric (due to Mac Fadden (TH [Ml [15] ) and the essentially imitative 
behavior among these agents (due to Brock and Durlauf [T71 [T5] ). 

Even though fundamental experiments dealing with social networks may constellate modern society 
analysis (i.e. the paradigmatic Milgram experiment of the sixties [35]), a real breakthrough in our 
understanding of network structure inside modern societies has been achieved when Granovetter reversed 
the Chicago school of social-psychology showing how a person with built weak ties (which were previously 
seen, at individual level, as ancestors of depressive states) was much more able to adapt its behavior to the 
social fitness due to the much broader amount of available information: in particular he noticed that these 
weak ties may often carry information of little significance (but not redundant as in highly clusterizcd 
community of similar agents linked by strong ties) , however they allow a primarily transmission of new 
information across otherwise disconnected cluster of the social network (allowing great potential benefit 
by these bridges). 

Two decades after this achievement. Watts and Strogatz, trough a mathematical technique (rewiring) 
have been able to display the Milgram results (sometimes known as "six degrees of separation" ) by which 
they understood that social (as well as others, i.e. biological [I][S]) networks can not be described by 
purely ordered or purely random graphs (i.e. Erdos-Renyi ones (5] [36]) due to correlations among nodes 
which allow for a much faster transmission of information (real graphs show high degree of cliqueness 
[371110]). pioneering a quantitative approach to these new networks, nowadays often called "small worlds" . 

In a different but related context, Mac Fadden has shown how to infer a model for econometric estimation 
of binary decision making (reflecting accurately several real cases in social structures f31| ) by introducing 
fundamental dichotomic degrees of freedom inside each agent mirroring its personal attitudes (i.e. a bit 
string ^ of K entries fi ~ 1, A' where each entry represent an attribute, i.e. /it = 1 accounting for 
smoking such that ~^ = +1 states that the i agent smokes while ""'^ = states that he does not, 
and so on). Once indexed individuals by i, i ~ 1...N, and assigned an Ising spin to each individual's 
choice (Ti = +1 for the agreement or ai = —1 for disagreement, he chooses to exploit data by assuming 
a single particle model into a suitable external field hi (the "field" influencing the choice of i) which is 
a function of the vector of attributes . Since for the sake of simplicity attributes are taken as binary 
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variables, the whole theory can be described in terms of an effective one-body Hamiltonan H as 

N . N K 

i i fj. 

where /i is a scalar parameter ruling the overall intensity of the external stimulus (whose capabilities of 
influencing a generic agent i are encoded into its bit-string ^h)- This parametrization of hi correspond 
to what economists call a discrete choice model |34| , and shows a remarkable link between econometrics 
and statistical mechanics {HN(a;£^) can be seen as a random field Ising model): In fact discrete choice 
theory has the same variational flavour of thermodynamics as states that, when making a choice, each 
person weights out various factors such as his own gender, age, income, etc, as to maximize in probability 
the benefit arising from his/her decision. 

Despite this result, there exist many examples from economics and sociology where it has been observed 
how the global behavior of large groups of people can change in an abrupt manner as a consequence of 
slight variations in the social structure (such as, for instance, a change in the pronunciation of a language 
due to a little immigration rate, or as a substantial decrease in crime rates due to seemingly minor action 
taken by the authorities [IHIIII])- From a statistical mechanical point of view, these abrupt transitions 
should be considered as phase transitions caused by the interaction between individuals that can not be 
accounted by a pure one-body theory. Indeed, Brock and Durlauf have shown [12] how discrete choice 
can be extended to the case where a global mean-field interaction is present (providing an interesting 
mapping to the Curie- Weiss theory (CW) [UlEI), thus further highlighting the close relation existing 
between the econometric and the statistical mechanical approaches to these problems. 

Instead of introducing Brock-Durlauf approach (which is a systematic translation of the CW scenario in 
social sciences) we go one step forward following the subsequent generalization obtained by diving the 
ensemble of the N decision makers into clusters, due to Contucci and coworkers [El [20]: Introducing a 
general two-body Hamiltonian H^^a; J) as 

N N 

Hn{(j;J) = - Jimm - '^h^cTi, (1.1) 

i,l=l 1=1 

they went over by defining a suitable parametrization for the interaction coefficients Ju . Since each agent 
is characterized by k binary socio-economic attributes, the population can be naturally partitioned into 
2'^ subgroups, which for convenience are taken of equal size: this leads to consider a mean field kind of 
interaction, where coefficients Ju depend explicitly on such a partition as follows 

Ju = ijir^Sgg', if i e 5 and Z e 5', 

which in turn allows us to rewrite as 

NJ ■ 

i7jv(CT; J) = --^( E ^aa'^gmg' + E ^9^9) 

S,9' = l 9=1 

where nig is the average opinion of group g, namely rUg = Yl!l^(g-i)N/2^^+i 

This idea of partitioning the interaction matrix into clusters of similar agents can be extended to a natural 
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limit (that we work out here) such that the size of these clusters approaches zero in the thermodynamic 
limit (so to preserve each individual identity and uniqueness for each agent, or a measure of attitude 
fluctuations inside the original idea of clusters): interestingly this leads to an interaction matrix so that 

-^.-^^ee;, (1.2) 
/J 

and naturally collapses the concept of magnetization in spin glasses 1321 to the one of retrieval in neural 
networks g], ultimately switching frustration into dilution (as ^ £ [0, +1] instead of ^ € [—1, +1]). 



2 The model and its topology: graph theory 



In this section wc look at the population as a graph and wc study its topological properties: each agent 
i is represented by a node; couples of agents displaying a positive coupling Jij > are said to 

be in contact or to interact with each other and this is envisaged by means of a link between i and 
j, whose weight is just Jy (eft. eq |1.2p . This picture mirrors the idea that socio-economic relations 
between individuals or firms are embedded and organized in actual social networks, which follows from 
the seminal work by Granovetter. 



As anticipated, each agent i is characterized by a binary string ^i, which might be thought of as the 
codification of agent's attitude, either positive (1) or negative (0), towards a given issue. All strings are 
taken of length K and each entry is extracted randomly according to 

^(?r = +i) = ^, ^(cf = o) = ^, (2.1) 

in such a way that, by tuning the parameter a G [—1, +1], the concentration of non null-entries for the 
i-th string pi = X^^Cf can be varied; consistently, also the topology generated by the rule in Eq. 11.21 
is varied. In particular, when a — )■ — 1 the system is completely disconnected (and only discrete choice 
survives), while when a — > +1 each link is present, being = 1 for any couple (so to retain the 
fully Brock-Durlauf approach). As we will see, small values of a give rise to highly correlated, diluted 
networks, while, as a gets larger the network gets more and more connected and correlation among links 
vanishes. In agreement with our modelization intent, repetitions among strings arc allowed. 



The main topological features of the emergent network have been investigated in [1] [7] , where it was 
shown that the average link probability among two generic nodes is 

and that, for large enough N and K , with K growing slower than linearly with TV, the degree distribution 
is multimodal. Therefore, the average degree for a generic node reads as z = pN . 

Apart from these global, long-scale features, the model also displays interesting properties concerning 
correlation among links, as we are going to deepen. 
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Small-world (SW) networks are characterized by two basic properties, that are a large clustering coeffi- 
cient, i.e. they display subnetworks where almost any two nodes within them are connected, and a small 
diameter, i.e. the mean-shortest path length among two nodes grows logarithmical with N. While the 
latter requirement is a common property of random graphs |361 [5] , and is therefore satisfied also by the 
graph under study, the clustering coefficient deserves more attention. Several attempts in the past have 
been made in order to define network models able to display such a feature [36l O |38] . For instance, in 
their seminal work. Watts and Strogatz |38| introduced a rewiring procedure on links, which can yield 
the desired degree of correlation. As we are going to show, in our approach SW effects emerge naturally 
from the definition of patterns and from the rule in Eq. II. 2[ that is to say, interactions based on sharing 
of interests (i.e. non-null entries) intrinsically generate a clustered society. 



Before proceeding, we notice that the same property can be addressed in different ways: we can say that 
the graph exhibits a large transitivity, meaning that if i is connected to both j and k, then j and k are 
likely to be connected; in modern network theory we say that the graph displays a large "cliquishness" 
and we measure it by means of the so-called clustering coefficient [36l HO] ■ The local clustering coefficient 
for node i is defined as 

c. = -^^y (2.3) 

Z^[Z^ - 1) 

where Ei is the number of actual links present within the neighborhood of i, whose upper bound is just 
Zi(zi — l)/2, namely is the number of connections for a fully connected group of Zi neighbors nodes. 
Then, the (global) clustering coefficient for the whole graph simply follows as the average of Ci over all 
nodes. Usually, for a graph with average degree equal to z, having a large clustering means that c is 
larger than z/N = p, which represents the clustering coefficient c^^ for an ER random graph with an 
analogous degree of dilution. 

As for the graph under study, we found that [7] c~p+l/p— l/{z — 1) > p, where in the last equality 
we used p < z — 1, which holds when N is large enough and the graph topology non-trivial [7]. More 
generally, one can notice that for this graph, the Zi neighbors of i are all nodes displaying at least one 
non-null entry corresponding to any non-null entries of ^f, this condition biases the distribution of strings 
relevant to neighboring nodes, so that they are more likely to be connected with each other. Indeed, 
when Pi = 1, it is easy to sec that Ci — 1, to be compared with cf^, namely the average link probability 
for node i, which turns out to be (1 + a)/2 < 1; analogous arguments apply also for larger values of p 
[7] . A numerical corroboration can be found in Fig. 12. 1[ which shows that c > c^^ in a wide region of 
values of a corresponding to non-trivial networks, i.e. for a larger than the percolation threshold and 
smaller than the fully-connected threshold. The clustering effect is especially manifest in the region of 
high dilution, where, for graphs analyzed here, c is even two orders of magnitude larger than c^^. Of 
course, when a approaches 1, the graph gets fully connected and c — )■ c^^ — >■ 1. These results are robust 
as N is varied. 



Finally, we mention another quantity used in ecology and epidemiology to quantify the existence of 
correlations among links, that is the so-called assortativity coefficient [321 a network is said to 
show "assortative mixing" ( "dissortative mixing" ) on their degrees whenever high-degree vertices prefer 
to attach to other high-degree (low-degree) vertices. While assortativity is typical of social networks, 
dissortativity is often found in technological and biological networks. 

The assortativity coefficient r can be defined as a Pearson coefficient to measure the correlation between 
the coordination numbers at either ends of a link; the ER graph corresponds to r = [351 ■ The measures 
performed on the graph under study suggest a dissortative behavior (r < 0), which is corroborated by 
the quantity {z)z', representing the average degree over the nearest- neighbors of a node with degree z' , 
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Figure 2.1: Main figure: Ratio between the clustering coefficient c measured for the correlated network 
under consideration and the coefficient c^^ corresponding to a purely random graph displaying an 
analogous density of links, as a function of the parameter a and for different values of K; the volume N 
is kept fixed and equal to 2000. Inset: behavior of the average degree z with respect to a for the same 
set of realizations corresponding to the main figure. 



namely: 

JV-l 

z,, = ^ zP(z;z'), (2.4) 

z=Q 

where P{z; z') is the conditional probability that a link stemming from a node with degree z' points to 
a node with degree z. As shown in Fig |2.21 a decreasing behavior of z^' is consistent with a dissortative 
mixing. The reason of this behavior is clear to see: while nodes corresponding to strings with large p can 
connect to most other nodes, nodes with small p have basically no chance to connect to other similar 
strings. This gets more evident for a small, as the concentration of such strings is larger. Interestingly, as 
highlighted in [37| , dissortativity has significant effects on the resilience (see next section) of the structure 
itself: dissortativcly mixed networks are less robust to the deletion of their vertices than assortatively 
mixed or neutral networks. 

As first remarked in |24| , real social networks are not only characterized by a small- world topology, which 
basically means large clustering and small diameter, but they also feature a peculiar coupling pattern. 
In fact, not only the neighbors of a given node i are likely to be connected, but they form communities 
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Figure 2.2: Neighbor assortativity for different system sizes (from top to bottom = 2000, = 
5000, N = 8000, depicted in different colors) and different values of K/N (for each choice of N, from 
top to bottom, K/N = 25, K/N = 100, K/N = 250). Notice that for K smah different modes can 
be distinguished. The two insets show in detail the case N = 2000, K — 80 (upper panel) and N — 
8000, K = 320 (lower panel). The value of a is fixed and equal to —0.6. 



such that intra-group links are expected to be stronger than inter-group links. In this way weaker ties 
work as bridges connecting communities strongly linked up. Interestingly, analogous properties arc found 
also for Q{N,K,a), in fact, it is intuitive to see that nodes displaying very similar strings are likely to 
be intensively connected with each other, hence forming a group, while, each of them, separately and 
according to the pertaining string, can be weakly connected with other nodes/groups. 
In order to deepen the role of weak tics, we perform two percolation processes where links are deleted 
either dcterministically or randomly: Given the pattern of couplings, in the former case we delete those 
with magnitude lower than a given threshold t £ [0, 1] meant as a tunable parameter; in the latter 
case we progressively delete nodes in a random way. We call 1 — / the fraction of links erased (in the 
former case / is a function of u) and we measure the size of the largest connected cluster (see also [3] for 
more details). Results are shown and compared in Fig. 12.31 Interestingly , when weak links are deleted 
first the graph starts to be disconnected at a value of / rather small, and this is a signature that such 
links work as bridge. On the other hand, strong ties are highly redundant [5]. Indeed, starting from 
a connected graph, the first nodes to get disconnected are those with small p, to fix ideas, those with 
p — 1; among these, the ones with the non-null entry in the same position were completely clustered in 
the original network. As the threshold t is increased, more and more nodes get disconnected; most of 
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Figure 2.3: Size of the largest eomponent versus the fraction of hnks in the system. Several size are 
depicted N = 2000, N = 3000, N = 4000, N = 6000, all corresponding to the same value a = 100 
and 7=1. Notice that when weak ties are deleted first, the way the size of the largest component (s) 
increases with / is smoother with respect to the random delation and also that (s) is smaller than 1 at 
relatively large values of /. 

them remain isolated (typically those with p < l), however, some non-trivial components survive. Such 
clusters arc made up of very similar strings with a relatively large number of non-null entries (p > l) 
and are therefore all closely connected. 

It is worth noting that such strongly clustered components emerge just in the "critical region" , namely 
where it is possible to detect nodes bridging two clusters and which play as "brokerage" between distinct 
group; this is a strategic position since it allows access to a more diverse set of ideas and information. 
The notions of homogeneity within groups and intermediacy between groups form the basis for the theory 
of " structural holes" introduced by Burt [M] . 

We finally comment on the resilience properties of the network under study, which can be as well inferred 
from the analysis on percolation processes. According to the situation, the stability of the network can 
be defined as its ability to remain connected or to still exhibit a giant component, under edge removal. 
In the former case, if weak links are the most prone to failure, our correlated network performs rather 
badly. Conversely, if we are interested in the maintenance of a macroscopic connected component, given 
that weak links arc the first to be deleted, our correlated network performs definitely better, as the 
percolation threshold grows slowly with N (see also [3]). 
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3 The model and its relaxation: stochastic dynamics 



We saw that the general structure of the Hamiltonian obeys 

N,N K N 

HN{a;0 = ^ E E^f^^'^J +E'**^- (3.1) 

i<j fl — l i 

It is worth recalhng that the attributes (^) are drawn randomly once for all, and so ar treated as 
quenched variables: this does not mean that a particular agent does not evolve in time changing his 
attribute distribution, but that, overall, one agent may switch to another and viceversa as far as the 
global attribute distribution is kept constant. 

The state of the system at this time is given by the average of all its building agents (such that we can 
introduce a "magnetization" as their average m = N~^^^ Ui), each of which evolving time-step by 
time-step via a suitable dynamics: 

Following standard disordered statistical mechanics approach |16j we introduce the latter accordingly to 

a.,{t + 1) = sign (tanh(/3</j,(i)) + m{t)) , (3.2) 
where (pi{t) is the overall stimulus felt by the i-th agent, given by 

N 



it)=N-'J2Jv'^3it) + Ht), (3.3) 



and the randomness is in the noise implemented via the random numbers rji, uniformly drawn over the 
set [—1,-1-1]. P rules the impact of this noise on the state ai{t + 1), such that for /3 = oo the process is 
completely deterministic while for /3 = completely random. 

In the sequential dynamics we are introducing, at each time step t, a single agent It -randomly chosen 
among the N- is updated, such that its evolution becomes 

P[ai, {t + 1)] = 1 (1 + ai^ it) tanh(/3^z, (t))) , (3.4) 

whose deterministic zero- noise limit is immediately recoverable by sending /3 — > oo. 
If we now look at the probability of the state at a given time t + 1, Pf+i((T), we get 



1^1 1^1 



where we introduced the N flip-operators Fi, i e (1, N), acting on a generic observable 0(cr), as 

Fi$(cri, -l-CTi, ...,crjv) = $((Ti,...,-CTi,...,CrAr) (3.6) 

such that we can write the evolution of the network as a Markov process 

Pt+iim) = J2w[m;m']ptim'), (3.7) 

rn' 

N 

N 



1 ^ 

W[m; m'] = 5rn,m' + T7 E (^K-P'i"l)(5„,^Fm' " w^{m)5, 
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with the transition rates Wi{m) = -^[1 — (7^ ta.nh{/3(pi)]. 

As the affinity matrix is symmetric, detailed balance ensures that there exists a stationary solution 
Poo{m) such that (restricting hi{t) ^ hi G M Vi G (1, A^)) 

W[m,m']Pao{'m') = W[m' ,m\Pao{'m). 

This key feature ensures equilibrium, which implies 

P N K N 
ij fj, i 

namely the Maxwell-Boltzmann distribution [TSIIS] for the Hamiltonian (|3.ip . 

In absence of external stimuli, and skipping here the question about the needed timescales for "thermal- 
ization" , the system reaches an equilibrium that it is possible to work out explicitly and that reproduces 
all the features stressed in the first section (as we are going to show). 

For this detailed balanced system furthermore, the sequential stochastic process (|3.2p reduces to Glauber 
dynamics such that the following simple expression for the transition rates Wi can be implemented 

W,(m) = (l + exp(/3Ai/(a,;0)) \ AH{a,-O^H{F,m-0-H{m-0- (3.9) 



4 The model and its equilibrium: statistical mechanics 

In the previous section we showed that, if the affinity matrix is symmetric (i.e. Jij = Jji), so that 
detailed balance holds, the stochastic evolution of our social model approaches the Maxwell-Boltzmann 
distribution (sec cq. (j3.8p ). which determines the thermodynamic equilibria. 

The latter are obtained by extremizing the free energy A{f3, a, h) = —Pf{f3, a, h) = u{f3, a, h)~f3~^s{(3, a, h) 
[u being the internal energy and s being the intensive entropy) that, as it is straightforward to check, 
corresponds to both maximizing entropy and minimizing energy (at the given level of noise /3, attribute's 
bias a and external influences h). Furthermore, and this is the key bridge with stochastic processes, there 
is a deep relation among statistical mechanics and their equilibrium measure Poo, hi fact 

Poo(a; e, h) (X exp(-/3if(a; t h)), A(/3, a, h) = -/3/(/3, a, h) = ^Elog^^ exp(-/?iJAr(a; t h)). 

a 

The operator E that averages over the quenched distribution of attributes ^ makes the theory not 
"sample-dependent": For sure each realization of the network will be different with respect to some 
other in its details, but we expect that, after sufficient long sampling, the averages and variances of 
observable become unaffected by the details of the quenched variables. 

Hence, once the microscopic interaction laws are encoded into the Hamiltonian, we can achieve a specific 
expression for the free energy, from which we can derive both the internal energy u(/3, a, K) as well as its 
related entropy s(/3, a, K): 

uiP,a,h) = -d0{PfiP,a,h))=N-\H{a;^,h)), (4.1) 
s(/3,a,/i) = fiP,a,h)+r'dpiPfil3,a,h)). (4.2) 
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The Boltzmann state is given by 

u^mcT,^) = ^ J E *(^;e)e-''^"''^'«\ (4.3) 

where the normalization Z is called "partition function" and the total average (<I>) is defined as 

($> =E[c^($(a,e))]. (4.4) 

We want to tackle the problem of solving the thermodynamics of the model trough the Hamilton- Jacobi 
technique p3] [6][9]pT]. 

Before outlining the strategy, some further definitions are in order here to lighten the notation (see also 
[7] for more details): taken g as a generic function of the quenched variables we have 



^ ^ /AT\ /T^'\ /I , „\ /I „\ N + K^lt~l, 



b 



where we summed over the probability P{1) that in the graph a number I of weighted links out of the 
possible N X K display a non-null coupling, i.e. £, ^ 0; this problem has been rewritten in terms of P{lb) 
and P{lc), where P{lb) is the probability that 4 (out of N random links) are active and analogously, 
mutatis mutandis, for P{lc) (on K random attributes): In fact, can be looked at as an x K matrix 
generated by the product of two given vectors like rj and Xi namely = rjiXfj,-, in such a way that 
the number of non-null entries in the overall matrix ^ is just given by the number of non-null entries 
displayed by rj times the number of non-null entries displayed by x- Hence, P{1) is the product of P{lb) 
and P{lc) conditional to Iblc = I. 

We can introduce now the following order parameters 



N 

M, 



1 ^ 

-Y^oji+i{a,), (4.6) 



N 

■I 

and the Boltzmann states w; arc defined by taking into account only I terms among the elements of the 
whole involved. 

Namely, Wi+i has only I + 1 terms of the type aa in the Maxwcll-Boltzmann exponential, all the others 
being zero: By these "partial Boltzmann states" we can define the average of the order parameters as 

JV-l 

(M) - ^(OM/- (4.7) 
/ 

We are now ready to show our strategy by defining the following interpolating free energy, depending by 
two interpolants, which can be though of as time and space in a mechanical analogy [23] [6] [9] [2T| 

^ ^ N K N N 

A{t, x) = -E log E exp E ^r^.''^^'^^ + E ^^""^ + ^ E ' (4.8) 

(7 ij i i 

where the random links have the same distributions of the standard as in any standard stochastic 
stability approach. Of course statistical mechanics is obtained when evaluating this trial free energy at 
t = — /3, X = 0. Let us work out the derivatives now: 

dA{t^x) ^_l l + a dM^^^l±a 

dt 2^ 2 ' ^ dx ^ 2 ' ^ ' 
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If we now introduce the following potential V(t, x): 

Vit,x) = i(i±^)2((A/2) - (M)2), (4.10) 

we can write the following Hamilton- Jacobi equation for the trial free energy 

dtAit, x) + i {d,A{t, a;)) ' + Vit, x) = 0. (4.11) 

When interesting at the replica symmetric regime (in a nutshell an approximation -which is widely 
believed to be correct even though not yet rigorously proved- in which we do not consider fluctuations 
of the order parameters in the large size of the population limit) we simply have to solve the free motion 
because replica symmetry means limTv-s-oo l^(i, x) = 0. The free field solution is given by the action in a 
generic point of the space-time plus the time-integral of the Lagrangian C = (ii^)^(M^)/2. Namely we 
can write ^ 

A{x,t) = A{xo,0) + I C{t')dt'. (4.12) 
Jo 

So we have 

A{xo,Q) = log 2 + (log cosh(xoC + /3/i.)) = log cosh(xo + /3/i) + ^-Y"- log cosh(/3/i) . (4.13) 

The time integral of the Lagrangian (as there is no potential) is simply i(i^)^(M^)t and the equation 
of motion is a straight line x{t) = xq + ^-^{M)t such that overall we can write a.t t = —/3 and x = — 

xo = ^(3 

A{P) = log2 + (i±^) logcosh (/3[(i±^(M) + h)]) + (i^) logcosh [(3h) - ^[l±^f{M)\ (4.14) 



Now we want to deepen the information encoded into equation (j4.14p : namely we want to recover by 
this solution all the theories of interaction introduced in section one in a quantitative way. 
Let us start forgetting the network, so with the two limits of MacFaddend independent particle model 
(a —1) and the pure Brock and Durlauf theory (a +1): 

lim ^(/3,a) = log2 + (logcosh(;3/i)), (4.15) 

a— ^ — 1 

hm ^(/3,a) = log2+(logcosh(/3AW + /i)))-^(A/2 ), (4.16) 
in perfect agreement with thermodynamics [51 118j . 

Note that when extremizing the free energy with respect to the order parameter (which is just one in 
both cases because in the former, as there is no network, the only decomposition trough cq. 14.71 is 
the independent sum of all the disconnected agents, while in the latter only one graph survives -the 
unweighted fully connected- and P{{M)) = 5{M — Mcw)i where with Mew we meant the standard 
CW magnetization), the response of the system is described by the hyperbolic tangent (nothing but the 
logit fit function used in econometrics): 

d(^M)A{P,a^ -l,h) = ^ TO = (tanh[/3/i]), (4.17) 
di^M)A{P,a^+l,h) = 0^ Mew ^ {t&nh[P {Mew + h)]). (4.18) 
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In all the other cases of interest (so for a ^ ±1) a distribution for weights on links is always present and 
weights are stronger for links among nodes that share higher amount of attribute similarity as shown in 
the graph theory analysis. 



Last step now should be achieving the critical line, i.e. by the control of the fluctuations of = 
•\/]V (M — M) . To fulfil this task it is straightforward to follow the approach of [U [7] (section four) , with 
the streaming now given by the transport derivative D = dt + {^^){M)dx [13]. Instead of performing 
these calculations which mirrors the ones detailed exposed in the paper [7] and depict a phase transition 
at /3c = ( ^^2°^ (^^ naively expected), we find instructive to bridge the two solutions (at /i = 0), 
namely the one obtained in section tree of [7] (eq. (4.25) in that paper) and (14.141) . 
Starting from the former, let us at first use the self consistency relation ( eq. 4.24 of the cited paper) to 
transform tanh"^ (^, P(OM/) = (i±ii)/3(M). This gives 

AiP, a) = log2 + (i±^) logcosh(/3(i±^)(M)) + ^^(Af^) - (i±^)^/3(M)2. 
The latter can be written exactly as 

A{l3,a) = log2 + (i±^)logcosh(/3(i±^)(M)) - ^^^(M^) 

by assuming (A/)^ = (Af^), which, in a nutshell, is sharply the request S* = (zero source limit) in the 
double stochastic stability approach and V ~ in the Hamilton- Jacobi approach. 



5 Summary and outlooks 



In this paper we tried to bridge over different aspects of modern quantitative sociology in a unifying 
perspective ultimately offered by a simple shift of the patterns in an Hopfield model of neural networks. 
The fundamental prescriptions of the Granovetter and Watts-Strogatz theories from topological view- 
point and Mac Fadden and Brock-Durlauf ones from social influences are found as different limits of 
this larger model, where, in proper (wide) regions of the parameters (/3,a, h), all these features can be 
retained contemporarily, offering a systemic view of social interaction. 

The idea that a model for the associative memory of the brain (though of as an ensemble of many inter- 
acting elementary agents) may work even for quantifying social behavior is in general agreement with 
the "universality" found in all these complex systems, however, while the neural networks share both 
positive and negative links (so to preserve a low synaptic activity, i.e. X]^^^ ~^ -^(0, 1)), this prop- 
erty is avoided in our context (for otherwise the role of weak tics could be played by highly conflicting 
peoples). As a consequence, despite the role of anti-imitative coupling is fundamental (as discussed for 
instance in [5]), it turns out that the greatest part of social interactions should be essentially imitative (as 
sociologists know well from long time). Moreover, by focusing on couplings generated from the sharing 
of common attributes, a small world structure naturally emerges, and, consistently with real networks, 
it is possible to detect strongly clustered sub-communities. 

So, our main goal when dealing with these techniques, is not discovering other hidden breakthroughs 
among pioneering speculations, but offering a quantitative, predictive model (and related methods for 
its solution) by which recover agreement with data and theories and improve society accordingly to our 
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will. 

In this sense, we think that now a great effort must be achieved in dealing with the inverse problem and 
its related data analysis, on which we plan to investigate soon. 
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